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We present a compact analytic formula for the two-loop six-particle maximally helicity violating 
remainder function (equivalently, the two-loop lightlike hexagon Wilson loop) in M = 4 supersym- 
metric Yang-Mills theory in terms of the classical polylogarithm functions Lit with cross-ratios of 
momentum twistor invariants as their arguments. In deriving our formula we rely on results from 
the theory of motives. 
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INTRODUCTION 

The past few years have witnessed revolutionary ad- 
vances in our understanding of the structure of scattering 
amplitudes, especially in A^ = 4 supersymmetric Yang- 
Mills theory (SYM). It is easy to argue that the seeds 
of modern progress were sown already in the 1980s with 
the discovery of the Parke- Taylor formula for the sim- 
plest nontrivial amplitudes: tree-level maximally helicity 
violating (MHV) gluon scattering. The mere existence 
of such a simple formula for a quantity which otherwise 
would have been prohibitively difficult to calculate us- 
ing traditional Feynman diagram methods signalled the 
tantalizing possibility that a great vista of unanticipated 
structure in scattering amplitudes awaited exploration. 

In contrast to the situation at tree level, it is fair to 
say that recent progress at loop level has mostly been 
evolutionary rather than revolutionary, driven primarily 
by faster computers, improved algorithms (both analytic 
and numeric), and software for multiloop calculations 
which has been made publicly available. Yet we hope 
that a great new vista of unexplored structure awaits us 
also at loop level in SYM theory. 

This paper is concerned with the planar two-loop six- 
particle MHV amplitude [l|, Q , which in a sense is the 
simplest nontrivial SYM loop amplitude. The known in- 
frared and collinear behavior of general amplitudeSjCon- 
veniently encapsulated in the ABDK/BDS ansatz 0,13]; 
determines the n-particle MHV amplitude at each loop 
order L > 2 up to an additive finite function of kinematic 
invariants called the remainder function R^\ Given the 
presumption of dual conformal invariance [H, 0] for SYM 
amplitudes (not yet proven, but supported by all avail- 
able evidence [3, 01 IH 0, IH ) , can depend on confor- 
mal cross-ratios only. Since there are no cross-ratios for 

(2) 

n = 4, 5, the first nontrivial remainder function is i? 6 . 

The same function is also believed to arise 

as the expectation value of the two- loop lightlike hexagon 
Wilson loop in SYM theory [l]l 14 1 (after appropriate 
subtraction of ultraviolet divergences, e.g. |lj|). Numer- 
ical agreement between the two remainder functions was 
established in [l], Q|. In a heroic effort, Del Duca, Duhr, 
and Smirnov (DDS) explicitly evaluated the appropriate 



Wilson loop diagrams to obtain an analytic expression 

(2) 

for Rq as a 17-page linear combination of generalized 



polylogarithm functions [la, LLTj (see also [18J). 

The motivation for the present work is the belief that 
if SYM theory is really as beautiful and rich as recent 
developments indicate, then there must exist a more en- 

(2) 

lightening way of expressing the remainder function i?g . 
Ideally, like the Parke- Taylor formula at tree level, the ex- 
pression should provide encouragement and guidance as 
we seek deeper understanding of SYM at loop level. 

(2) 

We present our new formula for i?g in the next sec- 
tion and then describe the algorithm by which it was 
obtained. 



THE REMAINDER FUNCTION R, 
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The remainder function Rq is usually presented as a 
function of the three dual conformal cross-ratios 

S12S45 S23S56 S34SQ1 

ui = , u 2 = , u 3 = , (1) 
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of the momentum invariants s 



though we will see shortly that cross-ratios of momen- 
tum twistor invariants are more natural variables. In 
terms of 



X: U^X , 



<1 + U 2 + U 3 - 1 ± y/A 
2uiu 2 u 3 



(2) 



where A = (u\ + u 2 + U3 — l) 2 — 4uiu 2 U3, we find 

3 / 1 

Re(ui,u 2 ,u 3 ) = ^2 i L 4( x t> x i) ~ 9 Li4 ^ ~ l l Ui ^ 

_i (y^LWi-iM)! +-J 4 + — J 2 + — . (3) 
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Here we use the functions 
1 



3 



8!! 



\og(x + x ) 



E ^^^g(x + x-r(h- m (x + )+h- m (x-)) (4) 
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and 



£ n (x) = - (Li„0) - (-l)"Li„(l/x)) 



as well as the quantity 



(5) 



(6) 



Note that in the Euclidean region where all Ui > 0, the 
xf never enter the lower half-plane and the x~ never 
enter the upper half-plane. The expression <j3j> is valid 
in the Euclidean region with the understanding that the 
branch cuts of Li„(a;^~) and Li„(l/a;~) are taken to lie 
below the real axis while the branch cuts of Li„ (x~ ) and 
\A n (l/x~l) are taken to lie above the real axis. (The 
quantities x\ x~ appearing as arguments of the logs are 
always positive.) In writing Q extreme care has neces- 
sarily been taken to ensure the proper analytic structure. 
For example one can easily check that J naively simpli- 
fies to ^log(x~ /x + ), but this relation only holds in the 
regions A > or u\ + + u 3 < 1. We caution the 
reader that any attempt to use any such naive relations, 
including the well-known relation between Li n (l/:r) and 
Li„ (x) , without careful consideration of the branch struc- 
ture, voids our warranty on (|3|). 

Besides its great simplicity, two notable features of (O 
which set it apart from the DDS formula are manifest 
symmetry under any permutation of the u i; and the fact 
that the expression is valid and readily evaluated for all 
positive Ui, in particular also outside the unit cube. 



DESCRIPTION OF THE ALGORITHM 

A Convenient Choice of Variables 

The DDS formula is expressed in terms of the classical 
polylogarithms Life as well as a collection of considerably 
more complicated multiparameter generalizations stud- 
ied by one of the authors [3] and defined recursively by 



G(a k ,a k -i,- . . ; z) 



G(a k -i,...;t)-^- (7) 

o i — afe 



with G(z) = 1, of which the harmonic polylogarithms 
familiar in the physics literature (20j are special cases. 

The parameters of the various transcendental functions 
which appear in the DDS formula involve not just the 
cross-ratios ([T]), but also the more complicated combi- 
nations 1 - m, (1 - Ui)/(1 - m - Uj), Ui + Uj, uf kl 



■ _ u k —ui±y/ (u k +ui) 2 -4iiju k ui 
jkl 2(1 — uj )u k ' 



1-Uj-U k +Ui±\/K 

This large collection of variables is redundant in an ineffi 
cient way, with many rather complicated algebraic iden- 
tities amongst them. 



Our computation is greatly facilitated by a judicious 
choice of variables which trivializes all of these algebraic 
relations. We choose to express the three itj by six vari- 
ables Zi valued in P 1 (with an SL(2, C) redundancy) via 



^23^56 Z1QZ34 
Ui = , U 2 = , U 3 



Z14Z36 



214^25 ' 



(8) 



where Zi 



Zi. One virtue of these coordinates is 



that A becomes a perfect square, so that the Uj k[ are 

rational functions of the Zij. (The v^ kl completely drop 
out as explained in the following subsection.) 

We anticipate that for general n the best variables for 
studying the remainder function will be the momentum 
twistors of [2l|. Indeed the z variables may be thought 
of as a particular simplification of momentum twistors 
which is valid for the special case n = 6 via the rela- 
tion (abed) cx ZabZacZadZbcZbdZcd- In terms of momentum 
twistors 



Ml = 



(1234) (4561) 
(1245) (3461)' 



(1456) (2356) 
(1256) (3456)' 



etc. (9) 



The Symbol of a Transcendental Function 

We define a function of transcendentality degree 
k as one which can be written as a linear combination 
(with rational coefficients) of fc-fold iterated integrals of 
the form 



T k = / dlogRi 



o dlogRk, 



(10) 



where a and b are rational numbers, Ri(t) are rational 
functions with rational coefficients and the iterated inte- 
grals are defined recursively by 



d log Ri o 
b / i-t 



o d\ogR n = 
dlogRi o-.-odlogRn-i) d\ogR n (t). (11) 



The integrals are taken along paths from a to b. When 
the Ri are rational functions in several variables the issue 
of local path independence (or homotopy invariance) is 
important (see [22|), and we have checked that i?g ' has 
this property. 

A useful quantity associated with is its symbol, an 
element of the fc-fold tensor product of the multiplicative 
group of rational functions modulo constants (see [22I 
sec. 3]). The symbol of the function shown in dTU|) is 



symbol^) = R x 



Rk, 



(12) 



and this definition is extended to all functions of degree 
k by linearity. 
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The group property for rational functions Ri modulo 
constants implies that 



R 



1 • • • 



(RaRb) 



■■■Rk=Ri---®R a < 
+Ri ■ ■ ■ (g) R b § 

i?i • • • <g) (cR a ) ® • • • R k = Rx ■ ■ ■ ® R a 8 



• • Rk 

■ Rk, 

■ Rk, 



(13) 
(14) 



for any constant c and rational functions Ri. 

The symbol of all functions appearing in the DDS ex- 
pression for the remainder function i?g are readily com- 
puted using the above definitions. For example the clas- 
sical polylogarithm functions 



Life (z) 



Li fe _i(t)dlogt, Lii(z) 



clearly have 

symbol(Li fe (z)) 



-(1-z) 



log(l - z) 

(15) 

(16) 



Finally, we note that knowledge of the symbol does 
not uniquely determine a transcendental function; there 
is an ambiguity consisting in adding lower degree func- 
tions multiplied by numerical constants of the appropri- 
ate degree. 



Constructing a Prototype 

Our first goal is to construct a simple function which 
has the same symbol S as that of the DDS formula. Inter- 
estingly we find that the v^ kl variables completely drop 
out of the symbol, even before switching to the z vari- 
ables. Although the full expression is too lengthy to re- 
produce here (see [13]); we nn d that S possesses a crucial 
and surprising symmetry property which in components 
(the subscripts here run over the 15 different z^'s) reads 



[3 a 



bed 



3 



bacd 



bde 



3 



bade] 



-(af)c,[)f)rf)=0. 

(17) 

It is known (see [24j]) that transcendental functions of 
degree less than four can always be expressed in terms 
of logarithms and the classical polylogarithm functions 
Life only. However an arbitrary function of degree four 
(or higher) does not have this property, so it is rather 
remarkable that the remainder function i?2 can be so 
expressed. 

The explanation for this miracle is Conjecture 1.19 
of [3], from which it follows that any symbol obey- 
ing (|17p can be obtained from a function involving only 
logarithms and the classical polylogarithm functions Life 
with k < 4. There is however no ab initio statement one 
can make about the arguments of those polylogarithms; 
they could in principle be arbitrary algebraic functions 
of the variables appearing in the symbol. 



The key to constructing the function from the knowl- 
edge of its symbol lies in identifying the symmetry prop- 
erties of the five kinds of terms which can appear, shown 
in the following table: 



Function 


A<E)A 




A<E)S 


SOS 








/ 


/ 


Li 3 (x)log(y) 






/ 


/ 


Li 2 (x) Li 2 (y) 


/ 


/ 


/ 


/ 


Li 2 (x) log(y) log(z) 




/ 


/ 


/ 


log(x) log(y) log(z) log(iu) 








/ 



Here A (S) stands for antisymmetric (symmetric), and 
a checkmark under A<g>A indicates that the function's 
symbol contains a piece which is antisymmetric under 
exchange of the first two arguments and antisymmet- 
ric under exchange of the last two arguments, etc. The 
function log(x) log(y) log(z) log(iu) has a fully symmetric 
symbol and hence only contributes to S®S. 

The above table suggests the procedure for construct- 
ing a function with the desired symbol S. First we isolate 
the A®A part of S and write down a linear combina- 
tion of terms of the form Li 2 (x) Li 2 (y) with that sym- 
bol, which is only possible if (|17[) is satisfied. After sub- 
tracting the symbol of this function from S we isolate 
the S ® A part which we fit with functions of the form 
Li 2 (a;) log(y) log(z). At the next step, A®S, we need both 
Li4(x) and L'i3(x) log(y), which can be isolated from each 
other because the former does not contribute to the part 
of the symbol which is antisymmetric under exchange of 
the middle two arguments, as evident in (|16l) . Ultimately 
we are left with terms which are completely symmetric 
under exchange of all arguments and can be fit via func- 
tions of the form log(x) log(y) log(z) \og(w). 

This procedure leads to a function of degree 4 whose 
symbol matches that of the DDS function using only 
classical polylogarithm functions. One interesting fea- 
ture which emerges is that the arguments of all polylog- 
arithms in the result are either "diagonal" cross-ratios 
such as (|8]l. which omit two diagonally opposite z's, or 
"edge" cross-ratios, which omit two neighboring z's, such 
as 



Z14Z23 

Z12Z34' 
Z14Z56 

^16^45 ' 



c 2 - 
,.+- 



gig z 25 

Z12Z56 ' 
Z25Z34 

Z23Z45 



Z36Z45 
Z34Z56 ' 
Z12Z36 

z\aZ23 ' 



(2) 



(18) 
(19) 



Fitting the Remainder Function R t 

(2) 

The prototype for R 6 is incomplete in two related 
aspects. First of all there is considerable ambiguity in 
the placement of branch cuts, which the symbol (since it 
involves only rational functions) is completely insensitive 
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to. For example one may be tempted to use the identity 
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1-aS 



= (x+/x-y 



(20) 



to replace J by | log(a; _ /a; + ) as mentioned above, but 
the two quantities have different branch structures in u- 
space. All such ambiguities are resolved by the physical 
requirement that the remainder function should be real- 
valued and smooth in the Euclidean region. 

Secondly it remains to determine the terms in ([3]) pro- 
portional to 7r 2 or 7r 4 , which do not contribute to the sym- 
bol. These terms can be fixed by several considerations 
such as checking the various Regge limits which were tab- 
ulated in [TtJ (see also 25 -29|), the requirement that i? 6 2 ' 



must vanish in the collinear limit R 6 (1 — u,u,0) = 0, 
or by numerically fitting to the DDS formula. We have 
also checked agreement with the leading term in the ex- 
pansion around the collinear limit computed in [30| . 



DISCUSSION 

We have presented in (J3]) a compact analytic expression 
for the two-loop six-particle MHV amplitude or Wilson 

(2) 

loop remainder function i?g in terms of the classical 
polylogarithm functions Life only. The simplicity of our 
result suggests that it is no longer outrageous to imagine 
the possibility of determining the remainder function for 
all n and any number L of loops, with an aim of matching 
onto the work of 30 34] at strong coupling. 

For n > 6 we imagine that the remainder function will 
continue to involve only the diagonal and edge momen- 
tum twistor cross-ratios, generalizing (|8]) and fp~8|) . It is 
also notable that our formula involves a particular com- 
bination (j4j of polylogarithm functions with special ana- 
lytic properties which is evidently very closely related to 
the Bloch-Wigner-Ramakrishnan-Zagier (BWRZ) func- 
tions (see for example (3||). These generalize very easily 
to any even degree 2L, suggesting a natural appearance 
in the L-loop remainder function. We anticipate that the 
consistency of collinear limits, and especially the system- 
atic expansion away from them recently studied in (30j . 
will tightly constrain the general form of i?i L '. 

Finally we cannot pass up the opportunity to point 
out that the Bloch-Wigner functions (the degree 2 case 
of the BWRZ functions) are known [35j to compute the 
volumes of hyperbolic tetrahedra, raising the fascinat- 
ing possibility of making a connection to recent work re- 
lating scattering amplitudes to volumes of polytopes in 
AdS 5 0,(33 (see also (HHH). 
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